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ABSTRACT Three types of nonlinear Schrodinger models with multiple 
length scales are considered. It is shown that the length-scale competition 
universally results into arising of new localized stationary states. Multista- 
bility phenomena with a controlled switching between stable states become 
possible. 
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1 Introduction 

The basic dynamics of deep water and plasma waves, light pulses in nonlin- 
ear optics and charge and energy transport in condensed matter and bio- 
physics 10, H y, Q is described by the fundamental nonlinear Schrodinger 
(NLS) equation 






(1.1) 



where 'ip{x,t) is the complex amplitude of quasi- monochromatic wave trains 
or the wave function of the carriers. The second term represents the dis- 
persion and L is the dispersion length (e.g. in the theory of charge (energy) 
transfer L^ — Ti /2m with m being an effective mass). The nonlinear term, 
y|'0PV'j describes a self-interaction of the quasiparticle caused either by 
its interaction with low-frequency excitations (phonons, plasmons, etc.) [^ 
or by the intensity dependent refractive index of the material (Kerr effect) 
[|6|. The function f{x) is a parametric perturbation which can be a local- 
ized impurity potential, a disorder potential, a periodic refractive index, 
an external electric field, etc. It is well known that as a result of competi- 
tion between dispersion and nonlinearity nonlinear waves with properties of 
particles, solitons, arise. One can also say that this competition leads to the 
appearance of the new length-scale: the width of the soliton C = L/\/V. 
The presence of the parametric perturbation f{x) introduces additional 
interplays between nonlinearity, dispersion and perturbations. In the re- 
cent paper by Bishop et al. H the concept of competing length-scales and 
time-scales was emphasized. In particular, Scharf and Bishop p, H| have 
discussed the effects of a periodic potential {f{x) — ecos{2Ti x/C,p)) on the 
soliton of the NLS equation, and shown on the basis of an averaged NLS 
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equation that for Cp/C^l ot >1 the periodic potential leads to a simple 
renormalization of the solitons and creates a 'dressing' of the soliton. But 
when Cp ~ C there is a crucial length-scale competition which leads to the 
destruction of the soliton. Another interesting example of the length-scale 



competition was provided by Ref. |10 where the authors showed that the 



propagation of intense soliton-like pulses in systems described by the one- 
dimensional NLS equation may be left practically unaffected by the disorder 
(when /(x) is a Gaussian (5-correlated process). This theoretical prediction 
has recently been confirmed experimentally using nonlinear surface waves 
on a superfluid helium film [[ll| . 

The goal of this paper is to extend the concept of the length-scale com- 
petition to the essentially non-integrable systems: to systems with nonlocal 
dispersion and to systems with unstable stationary states. 



2 Excitations in nonlinear Kronig-Penney models 

Wave propagation in nonlinear photonic band-gap materials and in periodic 
nonlinear dielectric superlattices [Q |3| consisting of alternating layers of 
two dielectrics: nonlinear and linear, is governed by the NLS equation 

+ w^S{x-Xn)\'ip{x,z,t)\'^'ip{x,z,t) = , (1.2) 

n 

where Xn — ni is the coordinate of the n-th nonlinear layer {£ is the distance 
between the adjacent nonlinear layers), and it is assumed that the width w 
of the nonlinear layer is small compared to the soliton width C within the 
layer. In this case the problem can be described by the nonlinear Kronig- 



Penney model given by Eq. (1.2). It was shown in Ref. uM that the field 
tjj{x,z,t) can be expressed in terms of the complex amplitudes ipnizjt) = 
ip{xm z, t) at the nonlinear layers. The complex amplitudes tpniz, t) can be 
found from the set of pseudo-differential equations 

C^K 2C^/t 



sinh ik tanh £k 

with periodic boundary conditions ipn+N = V'm where N is the number of 



layers. In Eq. (L3) the operator k is defined as kip = C~^y/—idt — C^ d'^ip. 
Equation (1.2) has an integral of motion - the norm (in nonlinear optics 
this quantity is often called the power) P = J_ \ip\'^dxdz. 

For the excitation pattern where the complex amplitudes are the same 
in all nonlinear layers, ipn{z, t) = ^(z, t), we get 



^-idt - e dl tanh ( ^ V-*9t - C' dl\ * - ^|*P* = . (1.4) 
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Equation ( |1.4[ ) clearly shows the existence and competition of two charac- 
teristic length-scales: the interlayer spacing i and the size of the soliton in 
the nonlinear layer C- When i? ^ C one can expand the hyperbolic tanhence 
and Eq. ( |1.4D takes the form of usual NLS equation 

(z9t + C'a,2)* + yl*l'* = 0. (1.5) 



In the opposite limit, when ^ ^ C and tanh ( -^\f^~'^ — 1, equation (1.4) 



takes the form 



V^^iS^^C^v]/ _ |vi/|2vi, = , (1.6) 



which reduces, for static distribution ipt'^ = 0), to the nonlinear Hilbert- 
NLS equation recently introduced in Ref. |l3|. It is noteworthy that in 
contrast to usual NLS solitons, the localized solutions of the nonlinear 
Hilbert-NLS equation have algebraic tails |]la . 

It is worth to note the close relation of the problem under consideration 
to the theory of the long internal gravity waves in a stratified fluid with a 
finite depth h (see e.g. [^) which are described by the equation 

atM + ia^^u + 2ua^u -f ra^u = , (1.7) 

where r(-) is the singular integral operator given by 

00 

—00 

{j).v. means the principal value integral). In the shallow water limit (h — > 0) 
the dynamics is described by the Korteweg-de-Vries equation, dtu+\d\ 



u~ 



2udxU = 0, while the Benjamin-Ono equation, dtu + Hd^u + 2udxU = 0, 
governs the water wave motion in the deep-water limit {h —> oo). Here 

oo 

{H f){x) — —p.v.J dyf{y)/{y— x) is the Hilbert transform. 

— oo 

Being interested in stationary states of the system we consider solutions 
of the form ^(z, i) = 0(z) exp(zAt), where A is the nonlinear frequency 
and 4>{z) is the real shape function. Since Eq. ( |l.3| ) is Galilean invariant, 
standing excitations can always be Galileo boosted to any velocity in z- 
direction. For the shape function ^(z) we obtain a nonlinear eigenvalue 
problem in the form 

VA~C^tanh(^^VA^C^) 0-^03^0. (1.9) 

Simple scaling arguments show that in the low- frequency limit (A^^/^^ — > 
0) the norm behaves in the same way as in the case of usual NLS equa- 



tion (1.5): P ~ \/A. When A£'^/(^ -^ oo the norm P is a monotonically 
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FIGURE 1. Power P of the stationary state ipn{z,t) = e' *</>(«) vs the nonlinear 
frequency A. Numerical results from Eq. (1.9) (dashed line), Pade approximation 
(full line) and the asymptotic relation P ~ A~^" as A — + oo (dotted line). 



decreasing function: P ^ 1/vA. From the analysis of Ref. [n4J follows 
that the norm P(A) is a nonmonotonic function with a local maximum at 
Am « 1.25^/^ (see Fig. |l|). Thus, the stationary states exist only in a finite 
interval, < P < P{Km) and for each value of norm in this interval there 
are two stationary states. This is an intrinsic property of the nonlinear 
Schrodinger superlattice system. 



Discussing the stability of the stationary states satisfying Eq. (l.£), there 
are two sources of instability to be considered: longitudinal and transversal 
perturbations. The perturbations of the first type are of the same symmetry 
with res pect to transversal degrees of freedom as the stationary states of 
Eq. (L9), while the second type of perturbations breaks this symmetry. It 
was shown in [|l4| that stationary states which correspond to the branch 
with dP/dK < are unstable due to the longitudinal perturbations. The 
states with A > (4^^/3^^) sin^(7r/A^) are, in their turn, unstable due to the 
transversal perturbations. Thus, one can expect stable stationary solutions 
for nonlinear frequencies satisfying the condition 



A < |^mm{-sm^(-),A„,}. 

In particular, this means that the stationary state 'ipn{z,t) 

neither exist in the case of only one nonlinear layer (i —>■ oo) nor in the 



„iAt, 



(1.10) 



[z] can 
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quasi-continuum limit {N —f oo). But in the latter case the system supports 
stationary states which are localized in both spatial directions (see Rcf. [Q 
for details). 



3 Discrete NLS models with Long-Range 
dispersive interactions 

In the main part of the previous studies of the discrete NLS models the dis- 
persive interaction was assumed to be short-ranged and a nearest-neighbor 
approximation was used. However, there exist physical situations that def- 
initely can not be described in the framework of this approximation. The 
DNA molecule contains charged groups, with long-range Coulomb inter- 
action (l/r) between them. The excitation transfer in molecular crystals 
[ [l7| and the vibron energy transport in biopolymers |]1^ are due to tran- 
sition dipole-dipole interaction with 1/r^ dependence on the distance, r. 
The nonlocal (long-range) dispersive interaction in these systems provides 
the existence of additional length-scale: the radius of the dispersive interac- 
tion. We will show that it leads to the bifurcative properties of the system 
due to both the competition between nonlinearity and dispersion, and the 
interplay of long-range interactions and lattice discreteness. 

In some approximation the equation of motion is the nonlocal discrete 
NLS equation of the form 

J-l-V'n + ^ Jn-m(V'm - V'n) + IV'nl^V'n = , (l-H) 

where the long-range dispersive coupling is taken to be either exponentially, 
Jn = Je~''l"l, or algebraically, J„ — J\n\~^, decreasing with the distance 
n between lattice sites. In both cases the constant J is normalized such 
that ^ J„ = 1 for all f3 or s. The parameters f3 and s are introduced to 

n=l 

cover different physical situations from the nearest-neighbor approximation 
(/? ^ oo, s ^ cxd) to the quadrupole-quadrupole (s = 5) and dipole- 
dipole (s = 3) interactions. The Hamiltonian H = ^ Jn-m\''Pn — V'mP — 

n,77i 

i ^ |-0„|^, which corresponds to the set of equations (|l . ll| ) , and the number 

n 

of excitations N = J2 l'0nP are conserved quantities. 

n 

We are interested in stationary solutions of Eq. ( |1 . 1 1[ ) of the form ipn (t) = 
(j)n exp(jAi) with a real shape function </>„ and a frequency A. This gives 
the governing equation for </>„ 
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FIGURE 2. N umber of excitations, A'^, versus frequency, A, found numerically 
from Eq. (1.12) for s = oo (full), 4 (dotted), 3 (short-dashed), 2.5 (long-dashed), 



2 (short-long-dashed), 1.9 (dashed-dotted). 



which is the Euler-Lagrange equation for the problem of minimizing H 
under the constraint N — constant. 

Figure shows the dependence A^ (A) obtained from direct numerical so- 
lution of Eq. (1.12) for algebraically decaying Jn-m- A monotonic function 
is obtained only for s > Scr- For 2 < s < Sc- the dependence becomes non- 
monotonic (of A/"- type) with a local maximum and a local minimum. These 
extrema coalesce at s = Scr — 3.03. For s < 2 the local maximum dis- 
appears. The dependence A^(A) obtained analytically using the variational 
approach is in a good qualitative agreement with the dependence obtained 
numerically (see fl^). Thus the main features of all discrete NLS models 
with dispersive interaction Jn-m decreasing faster than \n—m\^^'''' coincide 
qualitatively with the features obtained in the nearest-neighbor approxima- 
tion where only one stationary state exists for any number of excitations, 
N . However in the case of long-range nonlocal NLS equation (1.11), i.e. 
for 2 < s < Sen there exist for each N in the interval [A"i(s), A'^^u(s)] three 
stationary states with frequencies Ai(A^) < A2(A^) < A3(A^). In particular, 
this means that in the case of dipole-dipole interaction (s = 3) multiple 
solutions exist. It is noteworthy that similar results are also obtained for 
the dispersive interaction of the exponentially decaying form. In this case 
the bistability takes place for /3 < 1.67. According to the theorem which 
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FIGURE 3. Shows endpoints of the bistability interval for A'^ versus dispersion 
parameter s. For s — Scr the endpoints coalesce. Analytical dependence (full) 
gives Scr — 2.72. Numerical dependence (dashed) gives Scr — 3.03. 



was proven in [|l9|, the necessary and sufficient stability criterion for the 
stationary states is dN/dA > 0. Therefore, we can conclude that in the 
interval [Ni{s), Nu{s)] there are only two linearly stable stationary states: 
Ai(A^) and A3(A^). The intermediate state is unstable since dN/dA < at 
A = A2. 

At the end points (A(A^;) and A(A^„)) the stability condition is violated, 
since {dN/dA) s vanishes. Constructing the locus of the end points we obtain 
the curve that is presented in Fig. ^. This curve bounds the region of 
bistability. It is analogous to the critical curve in the van der Waals' theory 
of liquid- vapor phase transition [EOl. Thus in the present case we have a 
similar phase transition like behavior where two phases are the continuum 
states and the intrinsically localized states, respectively. The analog of the 
temperature is the dispersive parameter s(/3). 

The shapes of three stationary states in the interval of bistability differ 
significantly (see Fig. ^). The low frequency states are wide and continuum- 
like while the high frequency solutions represents intrinsically localized 
states with a width of a few lattice spacings. It can be obtained |Q that 
the inverse widths of these two stable states are ai ~ (iV/8J) = 

(7V/8J)'"^/^^"^'''*^ , as « ln(A^/J) with ^ = exp(l/s) being the char- 
acteristic length scale of the dispersive interaction which is defined as a 
distance (expressed in lattice spacings) at which the interaction decreases 
in two times. It is seen from these expressions that the existence of two so 
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FIGURE 4. Shapes of three stationary states for s — 2.5 and A'" = 3.1. The stable: 
A = 0.21 (full), A = 0.74 (long-dashed). The unstable: A = 0.57 (short-dashed). 



different soliton states for one value of the excitation number, N, is due to 
the presence of two different length scales in the system: the usual scale of 
the NLS model which is related to the competition between nonlinearity 
and dispersion (expressed in terms of the ratio N/ J ) and the range of the 
dispersive interaction ^. 

Having established the existence of bistable stationary states in the non- 
local discrete NLS system, a natural question that arises concerns the role 
of these states in the full dynamics of the model. In particular, it is of in- 
terest to investigate the possibility of switching between the stable states 
under the influence of external perturbations, and to clear up what type of 
perturbations can be used to control the switching. Switching of this type 
is important for example in the description of nonlinear transport and stor- 
age of energy in biomolecules like the DNA, since a mobile continuum-like 
excitation can provide action at distance while the switching to a discrete, 
pinned state can facilitate the structural changes of the DNA [|l) . As it was 
shown recently in ||22| , switching will occur if the system is perturbed in a 
way so that an internal, spatially localized and symmetrical mode ('breath- 
ing mode') of the stationary state is excited above a threshold value. 

We will in sequel mainly discuss the case when the matrix element of 
excitation transfer, Jn-mi decreases exponentially with the distance jn— m|. 
For /? = 1 the multistability occurs in the interval 3.23 < N < 3.78. It 
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FIGURE 5. Switching from continuum-like to discrete state for /3 = 1. The initial 
state (f>n has the frequency A ~ 0.31 and A'' — 3.6. The time evolution of |i/)„g (i)p 
when a phase torsion is applied to the center site with 6 — 0.261 (lower curve) 
and 9 — 0.262 (upper curve), respectively; inset shows time evolution of |i/'n(^)i 
for e = 0.262. 



is worth noticing, however, that the scenario of switching described below 
remains qualitatively unchanged for all values of /3 < 1.67, and also for 
the algebraically decaying dispersive coupling with 2 < s < 3.03. 

An illustration of how the presence of an internal breathing mode can af- 
fect the dynamics of a slightly perturbed stable stationary state is given in 
Figs. g| and |6|. To excite the breathing mode, we apply a spatially symmet- 
rical, localized perturbation, which we choose to conserve the number of 
excitations in order not to change the effective nonlinearity of the system. 
The simplest choice, which we have used in the simulations shown here, is 
to kick the central site ng of the system at t = by adding a parametric 
force term of the form 6 



^5{t)ipn{'t) to the left-hand-side of Eq. (1.11 



As can easily be shown, this perturbation affects only the site uq ai t — 0, 
and results in a 'twist' of the stationary state at this site with an angle 
0, i.e. V'no(O) — 4'no 6*^- The immediate consequence of this kick is, as can 
been deduced from the form of Eq. ( [l.ll[ ), that Jj (IV'noP) '^ill be positive 
(negative) when > (0 < 0). Thus, we choose > to obtain switching 
from the continuum-like state to the discrete state, while we choose ^ < 
investigating switching in the opposite direction. We find that in a large 
part of the multistability regime there is a well-defined threshold value 6th'- 
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FIGURE 6. Switching from discrete to continuum-like state for j3 = 1. The ini- 
tial state (j>n has the frequency A ~ 1.423 and N = 3.6. Same as Fig.H with 
9 — —0.490 (upper curve) and 8 — —0.491 (lower curve), respectively; inset 
shows time evolution of |'i/;„(t)| for 9 — —0.491 (only a part of a larger system 
is shown). 



when the initial phase torsion is smaller than 9th i periodic, slowly decaying 
'breather' oscillations around the initial state will occur, while for strong 
enough kicks (phase torsions larger than Ot^) the state switches into the 
other stable stationary state. 

It is worth remarking that the particular choice of perturbation is not 
important for the qualitative features of the switching, as long as there is 
a substantial overlap between the perturbation and the internal breathing 
mode. We believe also that the mechanism for switching described here can 
be applied for any multistable system where the instability is connected 
with a breathing mode. 



4 Stabilization of nonlinear excitations by disorder 



In this section we discuss disorder effects in NLS models. Usually the inves- 
tigations of disorder effects have been carried out on systems that are inte- 
grable - soliton bearing - in the absence of disorder. A common argument 
is that the equations, despite their exact integrability, provide a sufficient 
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description of the physical systems to display the essential behavior. How- 
ever, the more common physical situation is that integrability, and thus the 
exact soliton, is absent. A relevant example of such an equation is the two- 
dimensional (or higher-dimensional) NLS equation. The two-dimensional 
NLS equation is nonintegrable and possesses an unstable ground state so- 
lution which, in the presence of perturbations, either collapses or disperses 
(see e.g. H H). 

We consider a quadratic two-dimensional lattice with the lattice spacing 
equal to unity. The model is given by the Lagrangian 



L 



E h^ 



dt 



V'n- 



C.C. 



H 



(1.13) 



where 



H 



E(i^' 



n-\-1.7n 



i^n 



1^. 



n,m+l 



V^n/ 



^IV'n.r^ 



^n,m \V^n.7n \ 



(1.14) 



is the Hamiltonian of the system. In Eqs. ( 1.13D and (1.14) {n,m) is the 
lattice vector {n and m are integer). The first two terms in Eq. (1.14) cor- 
respond to the dispersive energy of the excitation, the third term describes 
a self-interaction of the excitation and the fourth term represents diagonal 
disorder in the lattice. Here the random functions en,m are assumed to have 
Gaussian distribution with the probability p(e„^m) — ^^ exp[— (e„^„i/?7)^] 

and have the autocorrelation function {tn^min' ,m') = V^^nn'^mm', where 
the brackets (....) denote averaging over all realizations of the disorder. 
From the Lagrangian (1.13) we obtain the equation of motion for the exci- 
tation function in the form 



at 

+ \i>m.nf4'rn,n + £m.n1pTn,n = 0- (1-15) 



Equation (1.15) conserves the norm N — Y^ IV'n.mP and the Hamiltonian 



H. 



We are interested in the stationary solutions of Eq. ( 1.15| ) of the form 

ipn,Tn{t) = 0„,„i exp(iAt) , (1-16) 

with a real shape function (j)n,m and a nonlinear frequency A. 

Equation ( 1.15| ) together with Eq. (1.16) constitute a nonlinear eigen- 
value problem which can be solved numerically using the techniques de- 
scribed in Ref. [^. The dependences N{A) in the absence and in the 
presence of disorder are shown in Fig. M. It has previously been shown 
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FIGURE 7. The norm A'^ versus nonlinear frequency A for various disorder 
strengths rj. Homogeneous case 77 = (solid line), rj — 0.04 (dotted line), 77 = 0.07 
(dashed-dotted line) and 77 = 0.1 (dashed line). 



[ p6| pT] , O, B8| that the linear stability of the stationary states in the dis- 
crete case is determined by the condition dN/dA > 0. Thus, in the case 
without disorder (solid curve in Fig. ^ the low-frequency (0 < A < A^^n = 
1.088) nonlinear excitations in the discrete two-dimensional NLS model are 
unstable. It is important that in the continuum limit (A —>■ 0) the norm 
A''(A) tends to the non-zero value Nc ~ 11.7. 

Other lines in Fig. show the dependence A^ on A for the stationary 
solutions of Eq. ( 1.15| ) in the presence of disorder. The results have been 
obtained as averages of 150 realizations of the disorder. Several new features 
arise as a consequence of the disorder. In the continuum limit (A —f 0) we 
no longer have N = Nc with dN/dA = 0. Instead we have A^ — > with 
dN/dA > signifying that the disorder stabilizes the excitations in the 
low- frequency limit. The disorder creates a stability window such that a 
bistability phenomenon emerges. Consequently there is an interval of the 
excitation norm in which two stable excitations with significantly different 
widths have the same norm. 

Furthermore, we see that the disorder creates a gap at small A in which 
no localized excitations can exist, and that the size of this gap apparently 
is increased as the variance of the disorder is increased. It is also clearly 
seen that as A increases (decreasing width) the effect of the disorder van- 
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ishes so that the very narrow excitations are in average unaffected by the 
disorder. It is important to stress that this is an average effect, because for 
each realization of the disorder the narrow excitation will be affected. The 
narrow excitation will experience a shift in the nonlinear frequency equal 
to the amplitude of the disorder at the position of the excitation. 




-50 50 -50 50 -50 50 



t=250 




FIGURE 8. Evolution of an initial excitation of the norm A'^ = 10.4402 without 
disorder (upper part) and with disorder strength r; = 0.1 (lower part). 



The bistability we observe in Fig. occurs due to the competition be- 
tween two different length scales of the system: one length scale being 
defined by the relation between the nonlinearity and the dispersion, while 
the other length scale being defined by the disorder. A similar effect was 
observed by Christiansen et al. |g9| for the one-dimensional discrete NLS 
equation with a quintic nonlinearity. The latter is quite natural because 
as it is well known (see e.g. [30 ) the properties of the two-dimensional 



NLS model with a cubic nonlinearity are similar to the properties the one- 
dimensional NLS equation with a quintic nonlinearity. 

Having studied the stationary problem it is vital to compare the results 
to full dynamical simulations. Therefore we carry out a numerical exper- 
iment launching a pulse in a system governed by Eq. ( 1.15| ). Specifically, 
stationary solutions (1.16) of Eq. ( 1.15| ) with A = 0.14 (after reducing the 
amplitude of these solutions by 5%) were used as initial conditions of the 
dynamical simulations. Examples of the described experiment arc shown 
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in Fig. pi. As is seen the pulse behavior in the absence of disorder and in 
the presence of disorder (we presented here a reahzation corresponding to 
the disorder variance rj = 0.1) differs drastically. While the pulse rapidly 
disperses in the ideal system (the contour plot for t = 250 is absent because 
the pulse width is of the system size), the process is arrested in the disor- 
dered system. After some transient behavior the excitation stabilizes and 
attains an approximately stationary width. The dynamical simulations thus 
support the conclusion that otherwise unstable excitations are stabilized by 
the presence of disorder in the low frequency limit. 

Analytical theory of soliton states in disordered NLS models based on the 
collective coordinate approach and on the Rice's theorem from the theory 
of random processes ISl) is presented in Refs. 129] and [p2[. 



5 Summary 

In summary we have shown that the presence of competing length scales 
leads to multistability phenomena in nonlinear Schrodinger models. We 
have analyzed three types of the NLS models. The nonlinear Schrodinger- 
Kronig-Pcnney model presents an example where two competing length 
scales exist: the width C of the soliton in the nonlinear Schroodinger equa- 
tion and the interlayer spacing £. Due to the interplay between these two 
length scales the localized stationary states exist only in a finite interval of 
the excitation power. Two branches of stationary states exist but only the 
low-frequency branch is stable. 

In discrete nonlinear Schrodinger models with long-range dispersive in- 
teractions there exist three types of length scales: the soliton width, the 
lattice spacing and the radius of the dispersive interaction. Here the compe- 
tition of the length scales provides the existence of three branches of station- 
ary states. Two of them: low-frequency branch which contains continuum- 
like excitations and high-frequency one with intrinsically localized excita- 
tions, are stable. It is shown that a controlled switching between narrow, 
pinned states and broad, mobile states is possible. The particular choice of 
perturbation is not important for the qualitative features of the switching, 
as long as there is a substantial overlap between the perturbation and the 
internal breathing mode. The switching phenomenon could be important 
for controlling energy storage and transport in DNA molecules. 

Considering nonlinear excitations in two-dimensional discrete nonlinear 
Schrodinger models with disorder it was found that otherwise unstable 
continuum- like excitations can be stabilized by the presence of the disorder. 
For the very narrow excitations the disorder has no effect on the averaged 
behavior. Bistability which was observed in this case is very similar to the 
bistability that occurs in nonlocal NLS models. Here the bistability arises 
on similar grounds because of competition between the solitonic length 
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scale and the length scale defined by the disorder. 
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